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ABSTRACT 

We calculate the rate at which dark matter halos merge to form higher mass systems. 
Two complementary derivations using Press-Schechter theory are given, both of which 
result in the same equation for the formation rate. First, a derivation using the prop- 
erties of the Brownian random walks within the framework of Press-Schechter theory 
is presented. We then use Bayes' theorem to obtain the same result from the standard 
Press-Schechter mass function. The rate obtained is shown to be in good agreement 
with results from Monte-Carlo and N-body simulations. We illustrate the usefulness of 
this formula by calculating the expected cosmological evolution in the rate of star for- 
mation that is due to short-lived, merger-induced starbursts. The calculated evolution 
is well-matched to the observed evolution in ultraviolet luminosity density, in contrast 
to the lower rates of evolution that are derived from semi-analytic models that do 
not include a dominant contribution from starbursts. Hence we suggest that the bulk 
of the observed ultraviolet starlight at z > 1 arises from short-lived, merger-induced 
starbursts. Finally, we show that a simple merging-halo model can also account for 
the bulk of the observed evolution in the comoving quasar space density. 
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1 INTRODUCTION 

If the matter content of the universe is dominated by 
cold dark matter, then dark matter halos associated 
with g alaxi es are expected to form hierarchically. Bond 
et al. ( 1991 ) have shown that a mo re rig orous treatment 
of the work of Press & Schechter ( 1974 ) can be used to 



in the rate at which halos of some mass form. Conse- 
quently, there have been a number of attempts to es- 



obtain information about the build-up of structure. Cal- 
culations using Press-Schechter (PS) theory are based 
on a statistical analysis of the initial field of density per- 
turbations and do not include any non-linear dynami- 
cal effects, but nonetheless the results appear to agree 
well with dynamical, coll i sionless ("N-body") s imula- 
tions dLacey fc Cole 1994 ^omerville et al. 1998] ). 

The aim of this paper is to use PS theory to calcu- 
late the rate of formation of dark halos of specified mass. 
We compare the results with those of Monte-Carlo and 
N-body simulations, and we also use the calculations to 
argue that the bulk of the observed cosmological evolu- 
tion in star-formation and quasar activity is a reflection 
of the evolution in rate of dark halo formation. 

A standard application of PS theory is to calcu- 
late the mass function given a cosmic epoch of interest. 
In many cosmological situations we are also interested 



timate this rate using PS theory (Lacey & Cole 1993; 
Sasaki 1994). In this paper we use the theory of ran- 
dom walks within the PS framework to calculate di- 
rectly the halo formation rate for any cosmology. We 
also use Bayes' theorem to derive the same equation 
from the standard PS mass function. These complemen- 
tary methods provide important insights into how the 
formation of halos can be understood within PS theory. 

In section ^ we use a Monte-Carlo realisation of ac- 
tual Brownian random walks to calculate the formation 
epochs of the halos they represent. The distribution of 
these epochs is found to be in good agreement with the 
formula calculated previously. Following this we show 
that the formation time distribution of halos in a large 
N-body simulation (found using a standard friend-of- 
friends algorithm) is also well approximated by the PS 
result. 

We then discuss the relevance of this work to our 
understanding of the observed cosmological evolution 
in star-formation rate (SFR) and quasar space density. 
In particular we suppose that there is a causal link be- 
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tween the hierarchical formation of dark halos and the 
activation of both quasars and luminous bursts of star 
formation. The possibility of a connection between these 
two key evolving quantitie s has already been suggested 
by a number of authors (Shaver et al. 1996 ; Dunlop 
1997; poylc fc Terlevich 199S ). Mergers and intcrac- 



density of the filtered field at any point is given by a 
Brownian random walk with the variance of the filtered 



field being the ' time' axis ( [Peacock fc Heavens 1990 



tions have long been implicated in both luminous star- 
bursts and active galaxies (see the review by Barnes 
& Hernquist 1992 and refs therein) and merging galax- 
ies have been found in which both phenomena are ob- 
served ( Canalizo fc Stockton 1997 ; ^tockton, Canalize 



Bond et al. 1991| ). The mass associated with this point 
in space is given by the position of the first upcrossing of 
an absorbing barrier at 6 = Sdt). The probability den- 
sity function that a trajectory will have its first upcross- 
ing at a mass between M and M + dM, P{M\5c)dM 
is given by a solution of the diffusion equation with an 
absorbing boundary condition (Bond et al. 1991): 



& Close 1998; Brotherton et al 1999| ). Investigation of 
the galaxies found at rcdshifts around 3 by the Lyman- 
break method ( ^tcidcl et al. 1998a ) show that much of 
the ultraviolet light inferred to be due to star formation 
can be identified with individual galaxies with star for- 
mation rates of possibly up to 1000 Mq per year (Steidel 
et al. l |998bD . These galaxies are strongly clustered and 
hence are likely to be the progenitors of cluster galax- 
ies at the present epoch (Steidel et al. 1998a, b; Adel- 
berger et al. 1998): at high redshifts there is therefore 
direct evidence for a link between the high level of star 
formation and host galaxies which are inferred to be 
merging into higher-mass galaxy systems. The detailed 
physics of both quasar- and star-formation is compli- 
cated and not understood, and in this paper we discuss 
only the role that the cosmological variation in halo for- 
mation rate might have in determining the evolution in 
observed quantities such as the ultraviolet luminosity 
density arising from star formation. 



P{M\5c)dM 



6c 



2^M 



dalj 



dM 



dM{\) 



where (t\j is the variance of the filtered field which for 



a sharp k-space filter is given by: 



>M 



1 r^/i? 



(2) 



where P{k) is the power spectrum of the initial random 
field. We can use equation || to obtain the comoving 
number density of halos of mass M at a given time t, 

n{M\t): 



n{M\t)MdM 



2p 



= pP{M\Scit))dM 
Sc{t) f 5S? 



(27r 



11/2^2 



■ exp 



M 



2^M 



duM 



dM 



dM. (3) 



This is the standard PS mass function. It should be 
noted that when using other filters this result is not 



gashima fc Gouda 1996) 



correct (Bond et al. 1991; Jedamsik 1995; Yano, Na- 



2 PRESS-SCHECHTER THEORY 

We now describe briefly the principles of Press- 
Schechter (PS) theory and calculate the number den- 
sity of halos at a particular epoch - the standard PS 
formul fi; a result which has been shown to be in good 
agreement with N-body simulations (e.g. Efstathiou et 
al. 1988). 

Dark halos are assumed to form by the non-linear 
gravitational collapse of initial density perturbations. 
We assume that the initial perturbations form a homo- 



To calculate the numbers of halos and their rate 
of formation we need to normalise the power spectrum. 
We require the variance in the density field when filtered 
with a top-hat filter of radius 8/i~^Mpc, as, to match 



the values deduced from X-ray clusters (Eke, Cole & 
Frenk 1996| ) 



geneous, isotropic Gaussian random field. In PS theory 



such a field is smoothed by convolving with a filter func- 
tion whose size is related to the mass M of halo in which 
we are interested. The fractional overdensity at any lo- 
cation is assumed to grow linearly until a critical over- 
density, (5c, is reached, when that location is considered 
to have collapsed into a dark halo of mass M, provided 
that the critical overdensity is not exceeded when the 
field is filtered on a scale corresponding to a larger mass. 
In the analysis that follows, instead of viewing the field 
as growing with time, we consider the field to be fixed 
and the critical overdensity 6c to decrease with time. 

If we use a sharp k-space filter of radius R, 
W{y;R) ^ (sin(r/i?) - {r / R) cos{r / R)) / {r / Rf , the 



For a top-hat filter, we can relate the mass of a 
halo to the filter radius, M = A/SirpR^. We can also 
calculate the critical overdensity for collapse for a uni- 
form spherical region; for a flat Qa ~ universe this is 
given by 6c = ol(\ + z). where a ~ 1.68 ( punn fc Gott 
1972). For an open f^A — universe, 6c(z) is given by 
Lacey fc Cole (1993|), and for a flat SIa ^ universe 



by Eke, Cole fc Frenk (1996). For a sharp k-space filter, 
the relationship between the mass and fil ter size is less 
obvious. Following Lacey and Cole ( |l993| ), we integrate 
the filter function over all space to obtain its 'volume 
of influence', which gives M = i/Air^pR^ . However, if 
we use the critical overdensity applicable to the top-hat 
filter used to normalise the power spectrum the k-space 
filter predicts a different number density of halos. We 
therefore choose the critical overdensity associated with 
a sharp k-space filter to be that which predicts the same 
number density of halos for both filters at the mass used 
to normalise the power spectrum. 
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3 DERIVATION OF THE HALO 

FORMATION RATE FROM PS THEORY 

We now use PS theory to obtain more information 
about the build-up of structure, in particular to cal- 
culate the rate at which dark halos of a particular mass 
form. We cannot calculate this simply by taking the 
time derivative of equation y as, in PS theory, halos 
are not only being continually formed at any particu- 
lar mass but are also continually being lost into halos 
of higher mass. We present new work to calculate the 
required halo formation rate by examining properties 
of the Brownian random walks invoked in standard PS 
theory. We use the fact that the rate of halo formation 
is equivalent to a probability density function in time. 
The analysis of Brownian random walks is important in 
many fields of pure and applied science and similar re- 
sults to those derived below can also be found in many 
introductory books on stochastic processes (e.g. Karlin 
& Taylor 1975). 

Initially, suppose that we are not interested in how 
the halo formed so we are not interested in the shape 
of the trajectory up until the epoch of halo formation. 
We wish to calculate the distribution of cosmic epochs 
at which halos of a given mass M are first formed, 
P(t|M)dt. In PS theory, this probability density func- 
tion is related by a variable transformation to P{5c\(7\j), 
the probability that a trajectory has its first upcrossing 
between 5c and 5c + d5c at a given value of a^^, and it 
is a formula for this that we now derive using the tra- 
jectories approach. We let 5{a1,j) be the position of the 



walk at aff and define: 



c(^m] 



max{(5(crjyj/ 



■,0<alj, 



< <^li}, 



<5diff (o-m) = ^maxlo-Af ) ^ ^i'^M) 



(4) 
(5) 



By considering the reflection of trajectories about the 
line 5 ~ m beyond the first upcrossing of this line, we 
can see that there is a one-to-one correspondence be- 
tween trajectories for which 5{a\j) > 2m — x and tra- 
jectories for which 5-^sx{'^\i) > rn and 5{a'^j) < x where 
X < m. For a pictorial representation of this correspon- 
dence see Fig. Q. This means that: 

PiSmaA'^li) > m,d{alj) < x\(j\j) 



P{5{al,) > 2m 



1 — erf 



xlcrli) 



2m 



(2^if)^/' 



, m > 0, m > x. (6) 



We differentiate with respect to x and then with respect 
to m, changing the sign, to get the joint density function 
for (5max(c^/) lying in the interval m, m-\-dm, and 5{a1,j) 



lying in the interval x,x 

f{m,, x\a]^,j)dm.dx 



dx. This gives: 



2 2m 



TT {a"^ 



M 



)3/; 



■ exp 



'{2m -x)' 
2^M 



dmdx, (7) 




2m-x 



Figure 1. The solid line represents a random walk consisting of 
70 equally spaced steps of (t|^ in the interval < o-|^ < 1.4. The 
random walk was chosen from a set of such walks and is such 
that <5max(l-4) > m and 5(1.4) < x. The dotted line represents 
the walk obtained by reflecting the initial trajectory about the 
line S = m after the first upcrossing of this line. Note that it has 
the property that (5(1.4) > 2m — x. Similarly, any trajectory which 
is such that <5(1.4) > 2m — x can undergo a similar reflection to 
obtain a trajectory with 5max(1.4) > m and <5(1.4) < x. There is 
a one-to-one correspondence between walks with 5max(o'|^) > m 
and S{a-j^j) < x and walks with (5(ct^,) > 2m — x. 



< rn, X < m. 

To obtain the joint density function for ^max(o'^) lying 
in the interval m,m -\- dm and (5diff(i7^) lying in the 
interval y,y -\- dy we note that: 

P{SmaA<^M) < a,5diff{alj) < b\crlj) 

f{m, m — y)dydm. (8) 

/o Jo 

From this we can deduce that the desired joint density 
is g(m,y)=f(m,m-y): 

9{m,y\(Tli)dmdy 

'2 m + y f—(rn + yY^ 



"^> 0,2/>0. 



exp 



2^M 



dmdy, (9) 



We now calculate g{m,\y = 0, alj) using Bayes' theorem 
to transform from equation O to the conditional proba- 
bility required and then setting y = 0. This is equivalent 
to setting y = in equation B and renormalising so the 
probability density function integrates to 1. Integrating 
g{m,y\a\,i)dmdy we see that: 
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9{yWli)dy ■ 



■ exp 



2<^M 



and applying Bayes' theorem we find that 



dy,y>0, (10) 



r I 2 N ™ + y 
ff(w|2/, O-A/j = 2 ^'^P 



'M 



2< 



dm, (11) 



rn > 0. 



It just remains to set y — and change the variable 
from m to Sc which we now take to be the position of 
the first upcrossing: 



PiScK,)dd, 



a 



Sc 

-ly- exp 

M 



lal, 



dSc 



(12) 



This is the probability that a trajectory has its first 
upcrossing between Sc and S,, + dSc at a given value of 

By changing the variables from aj^j to mass M and 
from Sc to time t, we obtain the probability that a halo 
of mass M formed between times t and t + dt. This 
gives the dependence on cosmic epoch of the rate of 
halo formation: 



P{t\M)dt = ^exp 



'M 



2-L 



dSr 



dt 



dt. 



(13) 



In the special case of a flat J7a = cosmology for 
which Sc — a(l + z), this time evolution is given by the 
simple formula: 



dn 



(x(l + z)^-^exp(-/3(l + z)^) 



(14) 



where (3 — a^ /{2a\[) is a function of the mass of halo 
and the power spectrum. 

We note that because the trajectories are Brownian 
random walks, all walks which pass through a given 
point can be thought of as new walks starting from that 
point. It is thus possible to apply coordinate changes to 
equation ^ and obtain the conditional probability for 
the mass distribution of progenitors of halos (Bond et 
al. 1990). Using Bayes' theorem and this result it is also 
possible to calculate the distribution of final collapsed 
masses attained by halos given the mass at an earlier 
time ( |Lacey fc Cole 1993| ). 

We can apply a similar argument to the result de- 
rived above for the time distribution of formation events 
and obtain the distribution of times at which a halo 
formed, given that it formed part of a larger halo at 
a known later time. If we consider walks starting from 
(cr^^,,(5^), where a\j, < a\j and S'^ < Sc (i.e. M' > M 
and t' > t) we obtain the required conditional probabil- 
ity: 

P{5c 



2 SI 

^M^Sc,cr 



M')dSc 



Sc-S' 



cr. 



- cr„ 



■ exp 



'M "M' 

The behaviour 



jSc-s'^r 

2('^M - <^l 



dSc 



(15) 



of a trajectory, having passed 



through a point related to the formation of a halo, thus 



provides mass and formation time distributions for the 
progenitors of the halo. As noted above, the form of 
the trajectory continuing from such a formation point 
is independent of the position of that point, and con- 
sequently the mass distribution of progenitors imme- 
diately prior to the formation event is independent of 
the formation epoch. Turning this argument around, 
we see that the effect of placing constraints on the pro- 
genitors of the halo immediately prior to the forma- 
tion event doesn't affect the distribution of formation 
times. We thus have the important result that, even if 
we only consider formation events resulting from merg- 
ers with specified progenitors, such as halos resulting 
from a merger of two approximately equally sized ob- 
jects, the distribution of times at which these mergers 
occur is still that given by equation 13. It is for this 
reason that we consider the cosmic variation in halo 
formation rate also to be a measure of the cosmic vari- 
ation in the rate of mergers between dark halos. We 
should note that because standard PS theory does not 
account for halo sub-structure, the derivation presented 
here excludes possible mergers between low-mass halo 
sub-units forming part of a larger collapsed halo. 



4 AN ALTERNATIVE DERIVATION 

We will now show that P{t\M)dt as given by equa- 
tion [13 can be calculated from the conditional prob- 



ability density function P{M\t)dM ( [Bond et al. 199lD 
using Bayes' theorem. Firstly, we wish to calculate the 
prior for Sc which is equivalent to asking the question 
'Given no knowledge of ct^, what is the distribution 
of upcrossing points in Jc?'- We note that all trajec- 
tories must have an upcrossing point of any given line 
S = Sc. Now, in any two equally sized intervals, dSi and 
dS2, there must be equal probability of such a crossing 
existing because the walk does not alter its form at dif- 
ferent Sc (all walks which pass through a given point can 
be thought of as new walks starting from that point). 
Thus, given no a priori information about a\j, all values 
of Sc are equally likely and we should assume a uniform 
prior for 6c. In order to be mathematically rigorous, we 
must make Sc bounded. We see later that we can remove 
these bounds without affecting the result. So we have 
that: 



P{Sc)dSc = 



dSa 



<-S„ 



^niin -1 ^c Jl ^max 

otherwise. 



(16) 



Applying Bayes' theorem and using equation n^, we can 
find the joint probability of Sc and afj-. 

Pi(7M,Sc)daljdSc = -rj —J r 

^" exp (-JSr) daljdSc. (17) 



(2^)1/2(4^)3/2 



2< 
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We can integrate this equation over all 5c to obtain the 
probability density function for a\j: 



M 



(2^)1/2(^2^)1/2(5, 



exp 



2-L 



exp 



6^ 



(18) 



and again use Bayes' theorem to obtain the conditional 
probabihty P((5c|ff/): 

6'. 



P{ScWM)dSc = -^- exp - 



'M 



exp 



In the limit as 6u 



^<i 



^< 



exp 



2^M 



2 

max 



(19) 



and (5, 

5] 



P{Sc\cr^i)dSc = -o- exp - 



a 



M 



max 
2 



2^1/ 



oo this becomes: 



ddc- (20) 



This is the probability that given a particular value 



of ctL, the first upcrossing at this cr^^ has probability 



P{5c\o'\i) of being between 5c and 5c + dSc- This prob- 
ability is identical to that calculated directly from the 
Brownian random walks in section pi 



5 COMPARISON WITH RESULTS FROM 
MONTE-CARLO SIMULATION 

As a check on the validity of equation |lj we have per- 
formed Monte-Carlo simulations to calculate random 
trajectories and determine the distribution of forma- 
tion epochs. Approximately 10^ random walks were con- 
structed, each consisting of 2^^ uniform steps in afj be- 
tween and a'^j equivalent to a mass of 1 .3 x 10^^ Mq for 
a standard CDM power spectrum with shape parameter 
r = 0.25 normalised to as — 0.64. Of these walks, 10'' 
were recorded with first upcrossing points at the final 
value of afj . The distribution of these upcrossing points 
was recorded in uniform bins in time and is compared 
with the expected distribution given by equation 14 in 
Fig H Good agreement is demonstrated. 



6 COMPARISON WITH RESULTS FROM 
N-BODY SIMULATION 

We present results from a simulati on run using the Hy- 
dra N-body, hy drodynamics code ( Couchman, Thomas 
fc Pearce 1995 ). The simulation used 128"^ dark mat- 
ter particles in a flat universe with Om — 1, J^a = 0, 
h = 0.5, and standard CDM power spectrum with shape 
parameter T = 0.25. The power spectrum was nor- 
malised to o-g = 0.64. Groups of particles were found us- 
ing a standard friends of friends (FOF) algorithm which 
links two particles together if the mean overdensity of 
a grid of particles with the same separation is greater 



X 

o 



■ Q.LO 



O 



O 




1 



2 3 4 



redshift 



Figure 2. Comparison of Monte-Carlo results with Press- 
Schechter predictions of the halo formation rate at fixed final 
mass. Open symbols show Monte-Carlo results of the forma- 
tion time of a halo of mass 1.3 X 10^^ Mq for a standard CDM 
power spectrum with shape parameter F = 0.25 normalised to 
(Tg = 0.64. The solid curve shows the prediction of equation n4 
at this mass. The results arc normalised to give a peak formation 
rate of 1. 



than that predicted at the moment of virialisation for 
spherical top- hat collapse ( Peebles 1980 ). 

In order to calculate the rate at which halos form we 
analysed results from the N-body run at 104 different 
times, separated by equal intervals of time. By com- 
paring the FOF results at each output time to results 
from earlier times we were able to build up a picture of 
the hierarchical growth of structure. Halos with at least 
half of the component particles not observed in a halo 
of equal or higher mass at an earlier time were recorded 
as being formed in the time interval between this out- 
put time and the previous one. There is a problem in 
this analysis that we may miss formation events - a halo 
formed in the time interval of interest may have already 
merged into a halo of higher mass when we analysed the 
simulation. However we can estimate the potential er- 
ror caused by this effect by calculating the maximum 
amount of halos which could have been formed from 
the data on all new halos. The distribution of new ha- 
los of between 45 and 55 particles is shown in Fig. g. The 
symbols mark the average of the minimum and maxi- 
mum mass which could have been involved in mergers 
per unit time. The minimum is just the recorded mass 
in new halos of between 45 and 55 particles, the max- 
imum is the mass in all new halos of greater than 45 
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Figure 3. Comparison of N-body results with Press-Schechter 
predictions of the halo formation rate at fixed final mass. The 
solid symbols show the N-body results for groups of between 45 
and 55 particles, corresponding to a mass of approximately 1.3 X 
10^^ Mq, see text for details. The solid curve shows the prediction 
of equation |L4| at this mass with parameters as given in the text 
and normalised to the N-body values. 



particles which could have passed through the interval 
of interest. The error bars mark the positions of these 
maximum and minimum points. Note that these error 
bars only denote the error in measuring the mass forma- 
tion rate from the N-body simulation results - they do 
not include errors intrinsic to the N-body simulation. 

In order to count enough formation events we had 
to use a reasonably large interval in mass, correspond- 
ing to between 45 and 55 particles. However, the curve 
predicted from equation 



14 



is similar within this mass 
range, and the error caused by this effect will be small. 
In Fig. ^ we plot the prediction from equation p^ at a 
mass equivalent to 50 particles which shows remarkably 
good agreement with the results of the N-body simula- 
tion. Note that we have not renormalised any of the 
parameters of the PS models other than by the self- 
consistent method given in section 0. One caveat we 
should note, however, is that we have only been able to 
test the agreement between the PS result and N-body 
simulation at relatively high halo masses. 



7 COMPARISON WITH PREVIOUS 
RESULTS 

We consider that the formation of a halo occurs when 
all of the mass in the halo is assembled. Blain & Longair 



model 


Hjvf 


f^A 


r 


h 


eg 


OCDM 


0.3 





0.15 


0.5 


0.85 


ACDM 


0.3 


0.7 


0.15 


0.5 


0.91 


SCDM 


1 





0.5 


0.5 


0.60 


rCDM 


1 





0.25 


0.5 


0.60 



Table 1. Table of parameters used for each of the 4 cosmo- 
logical models chosen. Here h = ii^o/lOOkms"^ Mpc~^, Qm = 
^^ Pmo- ^A = a 1 ^'^'^ r is the CDM power spectrum shape 
parameter. 



( 1993| ) and Sasaki ( 1994 ) also used the same definition 
of 'formation', and tried to calculate the rate of halo 
formation from the standard PS mass function (equa- 
tion y) and its derivative. However, they both made 
assumptions beyond standard PS theory. Blain & Lon- 
gair ( 1993 ) calculated the formation rate empirically 
assuming a form for the distribution of mergers which 
occur. Sasaki (1994) obtained an equation for the for- 
mation rate assuming the destruction rate (from which 
they then calculate the formation rate) for a power-law 
spectrum has no characteristic mass scale. These as- 
sumptions lead to different forms for the formation rate 
from that derived above where we have not made any 
assumptions beyond those of standard PS theory. 

Lacey & Cole (1993) defined the formation time 
as the time when the largest progenitor of a halo first 
contains at least half the mass of the halo. Using our 
definition of formation, we see that the distribution of 
such times is equivalent to the formation time distribu- 
tion of progenitors of mass greater than M/2 given that 
they arc the first such progenitors for a particular halo. 
The formation time distribution of progenitors of mass 
M given by equation ^ is only dependent on knowing 
that a progenitor of mass M is formed at some epoch. 
For instance by setting M — il/'/2 in equation |l^ we 
can calculate the distribution of times at which pro- 
genitors of mass M'/2 form, given that a progenitor of 
mass M'/2 is formed at some time. However, we do not 
know that the first subclump has mass M'/2, only that 
it lies in the range M' /2 < M < M' - this is a con- 
sequence of PS theory containing mass jumps. Lacey 
& Cole (1993) give a counting argument (their section 
2.5.2) which converts from a probability density in mass 
to a distribution in the number of progenitors in order 
to set the condition that we are only interested in the 
first halo to contain at least half the mass of the final 
halo. Without such an argument it is difficult to see 
how to distinguish the first progenitor from subsequent 
ones and so the results of this paper are not directly 
comparable to those of Lacey & Cole. 
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Figure 4. The formation rate of galaxy lialos as calculated using 
equation h3 for the FCDM cosmological model. Curves are plotted 
for halo masses in the range 10^"'^ — 10^^'" M© and are normalised 
to give the same peak formation rate. 



8 THE PREDICTED HALO FORMATION 
RATE 

We can now estimate the rate at which mergers occur 
to create halos of any particular mass for a number 
of cosmologies. We investigate 4 different cosmologies, 
based on those used by the Virgo consortium ( Springe] 
et al. 199q; Jenkins et al. 1998) summarised in Table 111 

Fig. kl shows how the functions vary with halo mass 
for the rCDM model. Curves are plotted at 4 different 
values of halo mass, 10^°'^ — 10^^ M©. Here it is easy to 
see the hierarchical build up of structure - low mass ha- 
los predominantly form first, with their peak formation 
rate occurring at higher redshift than that for higher 
mass halos. 

For each cosmology chosen, for a particular mass of 
halo, the redshift at which the formation rate peaks is 
highly dependent on the trg normalisation chosen. High 
values of ag, lead to early formation and the peak in 
the formation rate is at higher redshift. However by 
choosing a higher halo mass we can recover the same 
shape of curve. This is easy to see in the SCDM and 
rCDM cosmologies, where equation |lj applies, and the 
curve shape is only dependent on the value of (3 which 
can be kept constant by varying the halo mass with 
(78 ■ Fig. 1^ shows the mass merging rate for each of the 
cosmological models, but with halo masses chosen such 
that the models approximately agree on the redshift 
of the peak of halo formation. Viewed in this way, the 
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Figure 5. Comparison of the 4 different cosmological models. 
Because of the differing normalisations of the models the forma- 
tion rates are plotted for different values of the halo mass, chosen 
to provide the closest match between the models. The masses 
are: OCDM, IQI^-^Mq; ACDM, lO^^.Ojyig. gcDM, IOI^-IMq; 
rCDM, 10^^'^ M0. The curves are normalised to the same peak 
formation rate. 



OCDM model produces less evolution than the other 
models, but the other three have rather similar amounts 
of evolution. 



9 THE MERGER-INDUCED STAR 
FORMATION RATE 

A complete description of the formation of stars in 
hierarchically-forming galaxies requires us to under- 
stand mergers between dark halos, the fate of the bary- 
onic matter in such mergers, and the processes of cool- 
ing, star formation and feedback and stellar evolution 
that all affect the observed luminosity of a galaxy and 
its time variation. Some recent attempts to explain 
specifically the strong evolution in star formation have 
used ad hoc models of the rate of generation of new 
galaxies as a source function, which have then been 
combined with models describing the remaining pro- 
cesses and tested against the observed data in sub-mm, 
infrared, optical and ultraviolet wavebands (Madau et 
al. 1996, 1998). An alternative method is to simulate 
the merger histories of galaxies based on the work of 
Lacey & Cole (1993, 1994), add heuristic recipes for 
the formation of stars and predict quantities such as 
the galaxy luminosity function or the cosmic evolution 
in star formation rate (Kauffmann, White & Guider 
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doni 1993; |Cole et al. 1994| ; |Baugh et al. 1998|) . The 
success of these latter models has been the broad gen- 
eral agreement between the observations and the mod- 
els, in that evolution in the predicted star formation 
rate is seen, with a maximum rate at redshifts about 
unity (Madau et al. 1996, 1998; Baugh et al. 1998). But 
when considered in more detail, it is clear that those 
models do not reproduce the amount of evolution that 
is seen. That amount appears to be at least a factor 
of about ten (Lilly et al. 1996; Connolly et al. 1997; 
Madau et al. 1996, 1998), and the sub-mm data suggest 



O - 



that factors of about 30 may be required ( Hughes et al 
1998|)7puch large evolution in star formation rate would 
be consistent with the depletion of neutral gas that is 
inferred from the evolution in intervening absorption 
towards distant quasars (Pei & Fall 1995). The mod- 
els cited above appear to predict evolution only by a 
factor about 5. Semi-analytic models incorporating an 
additional evolving component of star-formation asso- 
ciated with short-lived starbursts have now been pro- 
duced ( Guidcrdoni et al. 1998D in which cosmological 
evolution by a factor (1 -I- z)^ is assumed, and these 
models do appear to give a better representation of the 
observed SFR. 

Our aim in this section is to provide a complemen- 
tary analysis to the above work by considering only the 
role of the cosmological evolution in halo formation rate. 
Observed luminous starbursts at low redshifts are prob- 
ably rather young, and simulations indicate that bursts 
of star formation associated with mergers between ei- 
ther disk or bulge/halo systems with a variety of relative 
masses all have lifetimes which are short (< 10^ years) 
(Mihos & Hernquist 1994, 1996). We have already ar- 
gued that the rate of halo formation should have the 
same cosmological evolution as the rate of mergers be- 
tween dark halos, and we should therefore expect any 
star-formation which is associated with mergers to have 
cosmological evolution consistent with the calculated 
evolution in halo formation rate. Similarity between 
the observed and calculated evolution would indicate 
that dark halo mergers are an important factor in star 
formation at high redshifts and we would then predict 
that most of the star formation in high redshift galaxies 
should be short-lived. 

The inferred star formation rates shown in Fig. 
are derived primarily from the observed luminosity den- 
sity at a rest-frame wavelength of 280 nm. Mod els of an 
evolvi ng galaxy after a burst of star formation (Bruzual 
& Cha :lot 1993| ) show that the ultraviolet light decays 
approximately exponentially with a timescale of about 
0.6 Gyr and it is this timescale which therefore domi- 
nates the light curve of short-lived starbursts. Thus, to 
compare the cosmic evolution in formation rate with the 
observed evolution we convolve the rate of halo forma- 
tion as derived above with an exponential function with 
this decay timescale. Note that this simple approach is 
not designed to provide a complete description of the 
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Figure 6. The observed mean comoving volume-averaged star 
formation rate as determined from tlie Canada-France redshift 
survey (Lilly et al. 1996; Madau et al. 1996) (solid circles), 
optical HDF data corrected for dust extinction (Madau et al. 
1996; Pettini et al. 1997) (crosses) and Connolly et al. (1997) 
(solid stars), extinction-corrected Lyman break galaxies (Stei- 
dcl et al. 1998c) (solid squares), sub-mm data (Hughes et al. 
1998) (open circle) and the star formation rate inferred from 
Ha surveys at zero-redshift (Gallego et al. 1995) (solid trian- 
gle) and at redshift 1 (Glazebrook et al. 1998) (open triangle). 
A Salpeter IMF and flat Qm = 1 cosmology has been assumed. 
The curves are the predictions of the FCDM model for masses of 
10l°0,10l'''5,10"°, 10^1-5 M0 normalised to the local star for- 
mation rate (see Section H). 



physics of star formation in these galaxies: merely to 
test the relative importance of the cosmic variation in 
merger rate. 

A sample of such curves are plotted with recent 
star formation rate data in Fig. 0: to be consistent with 
this previously- published data we adopt the Qm — 1, 
flA = cosmology and we show a range of halo masses 
for the FCDM model. Since the fraction of halo mergers 
which result in an observable starburst is unknown, we 
treat the normalisation of the models as being a free 
parameter and normalise the curves to the zero-redshift 
data-point. 

We see excellent agreement between the data and 
the cosmic variation in halo formation rate for any of the 
dark halo masses considered at redshifts up to about 1. 
At higher redshifts the data are very uncertain, and will 
remain so until corrections for optical extinction in star- 
bursts are better understood, but we expect the relative 
contributions from a range of halo masses to be impor- 
tant in influencing the observed form of the evolution at 
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Figure 7. The evolution predicted by combining tiie formation 
rate of different mass halos (convolved with an expected lifetime 
of 0.6Gyr) weighted by a Gaussian in dn/dlogM (solid line). The 
Gaussian is centred on a halo mass of lO'^'^'^ Mq and has u = 1.5. 
The evolution expected from halo formation events resulting in 
halos only of mass lO^"'^ Mq is also plotted (dotted line). Both 
curves are normalised to the data at 2 = 0. Data are plotted as 
in Fig. |. 



high redshift. As an illustration, wc combine formation 
rates for halos of a range of final masses, weighting the 
contribution from each mass by a Gaussian distribution 
in dn/dlogM. The relative normalisation for each mass 
is determined at z = using equation |[ Fig. shows 
the resulting evolution compared with the star forma- 
tion data, for a Gaussian weighting function centred on 
10^°'^ M0 with a logarithmic dispersion a = 1.5. Com- 
bining the curves in this way flattens the high redshift 
evolution as mergers creating lower mass halos become 
increasingly important. As expected, the < 2 < 1 
evolution is the same as that shown in Fig 0, but there 
is now much better agreement with the high-redshift 
data. A thorough treatment of this issue is beyond the 
scope of this paper, and would require additional mod- 
elling such as may be obtained from the semi-analytic 
approach. The illustration presented here is sufficient to 
show that the observed high-redshift behaviour can be 
reproduced by considering the star-forming galaxies to 
be dominated by short-lived merger-induced starbursts, 
and that it is the cosmic variation in halo formation rate 
that drives the evolution in star formation. 



10 QUASAR DENSITY EVOLUTION 

The comoving space density of quasars has also de- 
creased by a factor of order 100 from a redshift of 2 



to the present day (Boyle ct al. 1990). The quasar co- 
moving space density appears to have a maximum at a 



redshift about 2-3 (Schmidt, Schneider & Gunn 1995 
Shaver et al. 1996). Fig. H shows a compilation of mea 



suremcnts of the comoving space density of quasars se- 
lected in one decade of luminosity, shown as a function 
of redshift. The compilation shows quasars selected in 
radio, optical and X-ray wavebands: the consistency in 
their evolution is a strong argument that the form of the 
observed evolution is intrinsic to the quasar population 
and not due to either observational selection effects or 
to intervening obscuration (see also Shaver et al. 1996). 
The luminosity range that has been chosen is that in 
which the greatest amount of cosmological evolution is 
seen in the range < 2: < 2, and this therefore pro- 
vides the most stringent test of the model we discuss 
here. Previous authors have chosen instead to plot the 
comoving quasar luminosity density integrated over a 
wide range in luminosity (e.g. Boyle & Terlevich 1998): 
this produces a lower amount of evolution, in fact more 
comparable to the inferred evolution in SFR. We adopt 
the more pessimistic approach here. 

Our knowledge of the physics of quasar activation 
is even worse than our knowledge of the activation of 
intense bursts of star formation. We follow on from 
work by Efstathiou & Rees ( 1988| ) and Haehnelt & Rees 
(1993) and assume that merger events activate quasars, 
although we are not concerned with the actual mecha- 
nism of the activation. As for the induced starbursts we 
wish to establish whether the cosmic evolution in halo 
formation rate can account for much of the evolution in 
quasar numbers. 

We have little knowledge of the lifetime of quasars, 
but we can see that their lifetime must be short if the 
model is to reproduce the observed amount of evolu- 
tion. The association of a few quasars with starbursts 
might be taken as evidence that their lifetimes are sim- 
ilar, and the ages of extended radio sources are also 
thought to be short (< 10* years). We have therefore 
plotted the formation rate convolved with the same ex- 
ponential lifetime as used to model starburst evolution 
with the quasar data in Fig. & The normalisation has 
again been treated as a free parameter. Recent obser- 
vations of luminous quasars suggest that they live in 
massive host galaxies with a narrow range of host lu- 
minosities (McLeod & Rieke 1994a,b, 1995; Taylor et 
al. 1996; Hooper, Impey & Foltz 1997) and by infer- 
ence a narrow range of host masses. Here we compare 
the model at a single halo mass to the data. In this 
case the OCDM model does not produce enough evo- 
lution for any mass of halo and does not fit the data 
well. Reasonable fits can be obtained for the remain- 
ing three models, although all of them have a deficit in 
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Figure 8. The measured comoving space density of quasars 
is plotted for flat-spectrum radio-loud quasars (Sliaver et al. 
1996) witti P2.TGHZ > lO^''"'' WattHz-^sr-i (stars), optically- 
selected UVX quasars(Goldschmidt & Miller 1998; Boyle et al. 
1990) with -25.4 > Mb > -27.9 (open triangles), optically- 
selected high-redshift quasars(Schneider, Schmidt &c Gunn 1994) 
with 10'^* < Piy—cWatt < 10^® (open circles), and a compos- 
ite of ROSAT soft X-ray selected quasars(Boyle, Wilkes & Elvis 
1998; Boyle et al. 1994; Read et al. 1999; Schmidt et al. 1998) 
with lO^'''^ < Po.s,-2kevWatt < lO^*-'' (filled circles). The cos- 
mological parameters of the ACDM model have been assumed: 
the shape of the evolution shows little variation between the three 
cosmologies parameterised in Table 1. As each waveband selects a 
different subset of the overall quasar population, and also covers 
a slightly different range of bolometric luminosity, the normalisa- 
tion of each of the samples has been adjusted such that they lie 
on a consistent curve passing through all the data and then nor- 
malised to the X-ray data. The curve plotted is the predictions of 
the ACDM model for a halo mass of 10^^'^ Mq with normalisation 
of the curve adjusted to match the data (see Section hol). 



the amount of evolution between z = 2 and z = of 
about a factor 4. We have chosen to plot the ACDM 
model as, with its high ag, halos of masses 10^^-^ M© 
peak at z ~ 2.5. The SCDM model also fits reason- 
ably with a high value for halo mass, but can probabl y 
be excluded on other grounds (|Gawiser fc Silk 1998 ). 



The rCDM model fits with a value for the halo mass of 
j^Qio.e ]y[^. g^ value which is rather lower than would be 
implied by the observation that luminous quasars exist 
in massive host galaxies. The halo mass would be in- 
creased to 10^ ^'^ M© if (Tg were increased to 0.9: such 
high power-spectrum normalisations are indicated by 
the COBE data ( [Gawiser fc Silk 19981 ). Finally, we note 
that increasing h from 0.5 to 0.7 decreases the inferred 
halo masses by a factor about 3. 



11 CONCLUSIONS 

We have derived a simple formula for the rate at 
which structure merges to form new halos within Press- 
Schechter theory. This formula has been tested and 
shown to be in good agreement with the results of nu- 
merical simulations. 

We have also considered whether the strong cosmo- 
logical evolution in star-formation is primarily driven 
by the cosmic variation in the rate of halo formation. 
In order to produce a complete model of the evolution of 
merger-induced star formation wc need to know much 
more physics than we have used in this paper. However, 
it is expected that the evolution in the rate at which 
dark matter halos merge to form higher mass halos will 
have a large influence on the evolution of the merger- 
induced SFR. Examining Fig. || we saw that, although it 
is expected that a range of halo masses are important in 
producing star-bursts, for halo masses 10^° — 10^^'^ M© 
the observed drop in the SFR from redshifts z ~ 1 
to z '^ is independent of mass and the evolution 
predicted from any linear combination of these curves 
shows strong evolution. The similarity of this evolution 
to that observed is remarkable and, given that quies- 
cent star formation predicted from semi-analytic models 
does not provide enough evolution in this redshift range 
dGuidcrdoni ct al. 199^ ), indicates that merger-induced 
starbursts are extremely important for star formation 
at z ~ 1 and are perhaps the principal mechanism lead- 
ing to the observed star formation at high redshifts. 
This model does not preclude the undoubted existence 
of quiescent star-formation as well: it merely states that 
because the starburst component evolves so strongly it 
dominates at high redshift. At higher redshifts a more 
physically-motivated model is needed to deduce the rel- 
ative contributions of a range of halo masses, but we 
have shown that a simple combination of such a range 
can produce evolution wh ich is i n good agreement with 
the data of Steidel et al. (|l998c|) . 

The evolution of the quasar population at optical 
absolute magnitudes Mb ~ —26 is larger than that of 
the star formation rate, although the evolution of the to- 



tal quasar luminosity density is not (Boyle & Terlcvich 



199S). We have shown that, provided that quasar hfe- 



timcs are finite but short (around 0.6 Gyr), then merg- 
ers to form the halos of massive galaxies can reproduce 
the principal features of the observed quasar evolution, 
consistent with the observation that luminous quasars 
consistently live in massive host galaxies. But whilst 
the curve predicts very well all the basic features of the 
observations, it fails by about a factor of '^ 4 to pro- 
duce enough evolution between redshifts of and 2. In 
this sense the model does not provide a full explanation 
of the observed evolution, but as the PS theory does 
not include any non-linear dynamical effects we con- 
sider even this agreement to be surprisingly good. Pre- 
vious work has suggested that the expected cosmic vari- 
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ation in galaxy velocity dispersion (Carlberg 1990 ) or 
in halo circular velocity ( Haehnclt fc Rees 1993 ) could 



contribute to the quasar evolution, and only a mod- 
est amount of additional evolution would be needed. In 
the latter case we should expect to see some cosmo- 
logical evolution in the host masses associated with a 
given quasar luminosity. Some additional component of 



quasar luminosity evol ution ( Boyle ct al. 1990 ; Gold 
Schmid t fc Miller 1998| ) would also have the desired ef 



feet. Quasar lifetimes significantly longer than about 
1 Gyr would smooth out the evolution in these models 
and would produce worse agreement with the data. 

From the above comparisons between observations 
and the calculated rate of dark halo formation we con- 
clude that hierarchical merging is highly important for 
the strong cosmological evolution that is observed in 
both star formation and quasar activity. Determination 
of the masses of halos involved in star formation at high 
redshifts (e.g. Pettini et al. 1998) and quasars would 
provide good constraints on allowed cosmological mod- 
els and their parameters, and estimation of the ages of 
the star-forming systems at high redshift would be a 
direct test of the proposition that short-lived starbursts 
account for much of the observed star formation. 
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